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I. INTRODUCTION 



The electron-positron annihilation process is one of the basic hard scattering processes where the composite 
structure of hadrons is investigated. Especially, the two hadron production in inclusive e"*" e~ annihilation 
(THPIA), i.e. the process in which two hadrons are detected in the final state, attracts a lot of attention 
due to the possibility for both the experimental and theoretical studies of several fragmentation functions, in 
particular the chiral-odd Collins functions Since available experimental data 2, 3, 4] correspond to 

substantially different values of Q^, taking into account the evolution of fragmentation (in particular, Collins) 
functions becomes an important theoretical task. In turn, the studies of evolution arc intimately correlated to 
the investigation of QCD factorization. However, it is well-known that the factorization theorems do not hold 
for an arbitrary kinematic regime. Indeed, if a transverse momentum of the produced hadron is of the same 
order of magnitude as the large photon virtuality, pt ^ Q, the corresponding hadronic tensor can be factorized 
and expressed in terms of the integrated fragmentation functions (see, for instance, a detailed analysis of e+e~ 
annihilation in &]). Even if pr is much smaller than the large photon virtuality, but on the other hand, is much 
bigger than the characteristic hadronic size, Kqcd, factorization of the e^e~ annihilation cross section can still 
be proven with methods similar to those used for semi-inclusive DIS or Drell-Yan processes 7, 8, 9]. However, 
the regime for which the transverse momentum of produced hadron is of the order of Age d faces a conceptual 
problem how to identify the hard subprocess. Therefore, factorization becomes vague and requires a special 
care 0. 

Factorization in THPIA, starting from the seminal paper was only studied in the collinear approximation 
(see also 11 1). Besides, the procedure of integration over the transverse virtual photon momentum (see [l^ ) 
was restricted to Born approximation. 

Since the transverse momentum is of crucial importance for the Collins function, we perform a detailed 
analysis of its role in THPIA. We use ideas of 13|, [ij] where it was shown that integration over the transverse 
momentum of a produced /Lt~'~/x~ pair or a hadron in the Drell-Yan process provides an effective propagator of a 
highly virtual photon and generates the hard subprocess structure. This approach was recently generalized [15 1 
to include the case of the (weighted) transverse momentum average of semi-inclusive deep-inelastic scattering 
(SIDIS) with transverse momentum dependent fragmentation functions. 

In this paper, we develop these ideas further and present a new method of factorization which may be applied 
for any two-current process. In particular, we demonstrate the application of the method in the case of e+e~ 
annihilation when two hadrons, belonging to different back-to-back jets, are produced. 

To demonstrate the method we are proposing, we consider, at the first stage, the simplest case of the spin- 
independent fc_L-integrated fragmentation functions. We analyze contributions up to leading order as corrections 
to the hard part to obtain the evolutions of the corresponding fragmentation functions. Also, we briefly outline 
a way to extend our approach for the study of the Collins effects. 



II. KINEMATICS 



We consider the process e+(Zi) + e {I2) H{Pi) + H{P2) + X{Px) where positron and electron carrying 
momenta h and I2, respectively, annihilate into a time- like photon with momentum q = h + h for which = 



2 



is large. This time-like photon creates then two outgoing hadrons with momenta Pi {i = 1, 2) (these hadrons 
belong to two distinct jets) and an undetected bunch of hadrons with total momentum Px- For such kind of a 
process, it is convenient to introduce two invariants 

which are analogous to the Bjorken variable. Note that the variables zi and Z2 are the energy fractions of the 
detected hadrons in the e+ e~ center of mass system. They obey the following constraint due to the energy 
conservation: 

^ < 1. (2) 
Moreover, due to the momentum conservation there is a stronger constraint on zi and Z2'- 

21 < 1, Z2< 1, (3) 

Within this region, these variables can vary independently. 

To perform the Sudakov decomposition, we choose two dimensionless light-cone basis vectors: 

n; = (1/V2, Ot, 1/V2), = (l/x/2, Ot, -l/\/2), n* -n^l. (4) 

In this paper, we choose the kinematics such that the photon and one of the hadron have purely longitudinal 
momenta while the other hadron has both longitudinal and perpendicular momenta: 



V2 ^ Z1QV2 



The leptonic momenta are 



Q[1-COS02] , , Q[1+COS02] H fQ ■ n a 

2V2 2V2 "^"-^^^^ ^i,^(^-sm^2,0 

Q[l + cos6>2] , Q[l-cos6i2] 

2^2 2V2 ' V 2 



^2m = TTlK 'V + ^^^^ "m + ^2m. = ( S"^^2, ) , (6) 



where 62 is the angle between P2 and /i. This frame is called the " _L "-frame (or the perpendicular frame) 
[lij . Below, we will omit terms of order /Q^. Such a kinematics has advantage for analysis of the 
experimental situation, where the momentum of one of the produced hadrons is measured. 

In order to ensure that the two hadrons are in different jets, we introduce two different variables [lo| : 



Z^^i^^zi, U = ^^. (7) 



2Pi-g Pi-P2 
— — = zi, U = — 

Pi-q 

However, using eq. (O, the difference between U and Z2 is of the order of 

-1 



U = Z2 



l + Vi^/{zlQ^) , (8) 



and can be neglected in the leading order approximation which we are considering in this paper. 
The perpendicular projection tensor is defined as usual 

.9^^ = QtJ^v - qtj,qv + T^f^, (9) 

where the two normalized vectors q and T are constructed as 

* = |. (10) 
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For simplicity, we consider the unpolarized case. The differential cross section of the corresponding e^e 
annihilation is given by 

In terms of leptonic and hadronic tensors, we have 
where the hadronic tensor W^jy is defined as 



^^^=T.j ^^^S^^\q-P,-P2-Px){0\M0)\Pi,P2,Px){Pi,P2,Px\M0)\0)- (13) 
We can rewrite the part of the phase space corresponding to the detected hadron with momentum Pi as 

■d'Pi±. (14) 



Pi dzi 2 ■ 



(27r)3 2Ei (27r)3 2zi 

Because the leptonic tensor is independent of the hadronic momentum Pi, it is useful to introduce the 
averaged hadronic tensor, W^t^ 



j d^Pi^W^^\ (15) 



Finally, using and p5|). the cross section can be presented as 
Note that this averaging procedure produces the hard subprocess in the Born approximation 



^2 „^ ' = 77^ — i^'^KJ- (16) 



III. FACTORIZATION PROCEDURE: BORN DIAGRAM 

Let us now discuss factorization of hadronic tensors corresponding to the hadron production in e^- 
annihilation at low pT- In this case, as pointed out before, a serious conceptual problem is known to be 
associated with the difhculty to identify (or, in other words, to separate out) the hard subprocess. In [3], 
it was shown that in the case of Drell-Yan process a suitable integration over the transverse momentum pro- 
vides effectively the separation of the hard subprocess. We extend this approach to the factorization of the 
e"*" e~-annihilation hadronic tensor including leading order as corrections. Furthermore, the reproduction of 
the well-known DGLAP evolution kernel for both the quark and anti-quark fragmentation functions can be 
considered as a proof of longitudinal factorization. 

First, we consider a simple Born diagram depicted on Fig. [1] (a). The corresponding hadronic tensor reads: 

W^, = y d^fcid^fe 5(4) (fci-hfc2- 9) tr [7, 9(^)7^0(^1)] +(1^2), (17) 

where the four-dimension (5-function, representing the momentum conservation at the quark-photon vertex, will 
be treated as the "hard" part. The non-perturbative quark and anti-quark correlation functions 9(/c2) and 
G)(fci) are given by 

60^(^2)- j (0|^a(O|P2,Px.)(P2,PxJ#(0)|0), 

0a^(fci)--/ ^e-^'=-''(0|^MO)|Pi,Px,)(Pi,PxJV'a(ry)|0), (18) 
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where the underhned indices are the Dirac indices. The summation and integration over the intermediate unde- 
tected states are imphed. Note that owing to the delta function in (flT)) the hadronic tensor cannot be expressed 
through the factorized A:j_-integrated fragmentation functions. To express the corresponding correlation func- 
tions through the fragmentation functions, we suggest a factorization scheme inspired by the Ellis- Furmanski- 
Petronzio (EFP) approach We first consider the formal identical transformation: 

dz' 

d ki ^ d*ki dz[ 5{Pi /k^ — z[) = d ki ^ d{ki ■ h* — l/z[), 

dz' 

d^k2 ^ d*k2 dz'2 S{P+/k+ - z'^) = d^k2 5{k2 ■ h - l/z'2), (19) 

where the two vectors 

n* = fi. = (20) 

^' Pi-n*' ^ P2-n ^ ' 

have been introduced. Variables z[ may be interpreted as partonic fractions of the corresponding momenta of 
produced hadrons. 

As mentioned above, we treat the four-dimensional (5-function as the hard part of the corresponding tensor. 
This assumption will be justified below. Our analysis is limited by the study of the leading twist contributions. 
It means that we keep only the first terms of the expansion: 



<5W (fci + k2-q) = <5(4) (^ + ^-^+ 0{k 



q 

^2 



The hadronic tensor is now rewritten as 



r-^\ . (21) 



where 



^ / w / (S^Kf - - (t) " ' '''' 

e(z^) =^ J d^k2 5{k2 ■ n - 1/z^) e(fc2), (23) 
Q{z[) =^ / d^ki 5{ki ■ fi* - l/z[) e(fci) (24) 



Since we study spin-independent fragmentation functions, we first have to project the correlation functions (|23p 
and (j24p onto the corresponding Lorentz vector structures. Starting from the quark correlator function, we 
write 

e(4)=^ Jtr[7„e(4)]7a (25) 

Using and the integral representation of the (5- function in ([^5]) . the vector projection of the quark correlation 
function (sec (|25p ) can be represented in terms of the spin-independent fragmentation function as 

1 /■^^2g,A./4 ( dtdCd^U5{Xn- -r)5{^^)5'^^\L) 



4 7 27r 

tr[7„(0|^(e+,r,el)l^2,Px.)(P2,PxJVi(0)|0)]7a = A- (26) 



Here, the fragmentation function D{z2) can also be written as 



D{z'2) = ^ / dr e^^^^«"/^^ tr[7+ (0|V(0, T, 6x)|P2, Px.) (P2, Pxj^(0)|0)] , (27) 
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provided that the minus co-ordinate component £ is equal to . Note that this fragmentation function 
obeys the momentum conservation sum rules [i, la, isj : 



J dzzD{z) = 1. 

/l, s 



In a similar manner, we project the anti-quark correlation function: 



1 



tr [7„ Q{z[)\ 7„ 



"27' 



tr[7„ (0|V;(0)|Pi, Pxi) (Pi, PxJ^(?7+, ry±)|0)] 7a 



5(.i) 



Pi 



(28) 



(29) 



Taking into account Eqs. 



- (|29p . we bring the hadronic tensor in the following form: 



dz'2 



5(Pf M - q-)6{P}/z'2-q+)6^^\Vi/z[)tr ^.P^l^^Pi 



(30) 



Now, inserting P2 and Pj^ defined via the kinematical variables Z2 and zi (see (O) and calculating the trace 
in ([50]) . we get: 



zi / dz[D{z[)5{zi- z{) 



Z2 



Z2 



dz. 



{z 



2_ 

/\2 



D{Z'2)5{Z2 - z'2) 



(31) 



From PT|) . one can see that, though the hadronic tensor W/ii'' has a formally factorized form, it does not 
have much sense because of the two-dimensional (5-function. To eliminate it, we have to integrate over the 
perpendicular momentum or, in other words, to go over to the averaged hadronic tensor (|15p . 



4g^, I d2p,^5(2)(p^) [,^D{z,) 



D{Z2)_ 
Z2 



= -4g^, h^(2i)] 



D{Z2) 



Z2 



(32) 



On the other hand, as we will discuss below, integration over d^ Pi j_ will generate the effective diagram depicted 
in Fig. [1] (b). Indeed, after integration, one can represent the one-dimensional 5-function 5{z2 — z'2) as the 
imaginary part of the effective propagator: 



5{Z2 - 4) 



1 



^m- 



Z2 - Z2 



[9-^2/4]^' 



(33) 



The latter will reduce the diagram in Fig. [l](a) to the diagram plotted in Fig. [T](b). Notice that the appearance 
of this hard effective propagator justifies a posteriori the suggested generalization of the EFF factorization 
scheme. Namely, the four-dimensional (5-function should be treated as the hard part, and the entering parton 
momenta should be replaced, at the leading twist level, by their longitudinal parts. 

Let us now discuss the above-mentioned procedure from the viewpoint of the approach of averaging given by 



d^Pix = 2zi / / d*PiJ(p2)J 



2Pi -g 



z\ 



(34) 



Indeed, integrating over d^ P\ with the four-dimensional (5- functions (see ([HOI'), one can observe that the one- 
dimensional (5-function which is responsible for the mass-shell condition goes to the (5-function that can be 
understood as the imaginary part of some effective "propagator" with large photon virtuality, i. e. 



(35) 
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Diagrammatically, it corresponds to the case when the Born diagram transforms to the diagram plotted in Fig. 
[1] (b). The dashed Hne impHes the effective propagator or the factorization hnk. So, the averaged hadronic 
tensor can now be written in terms of the factorization hnk, 



wit' - / d^Pi^ ^ 1^ [,,D{z,)] [z, j ^i5(4)3m-' 



(36) 



The spuriously asymmetric form of (|32p or p6p with respect to the inter changing of zi and Z2 emerges because 
of the integration over the momentum Pi of one of the produced hadrons. Restoring the flavour dependence 
omitted above and by inserting into (fTB|) . we get: 

da{e+e-) 3 tt a^^^ ^ ^^^^ Y.^lD^{z,) D^{z,l (37) 



dzidz2d cos 92 2 



a.a 



where the only remaining asymmetry is reflected in the polar angle of the detected hadron. After integration 
over f?2, the result becomes completely symmetric: 

dcos02 , , = ^ j:-lD^iz,)D^iz2). (38) 

d zi a zod cos Oo Q ^ 

a.a 

Concluding this section, we would like to stress that the factorization of the Born diagram can be implemented in 
a similar way as in the case of Drell-Yan process [l2'|. In this case, we transform the two-dimensional integration 
d^ Pi _L in the phase space (see ()14p ) to the two-dimensional integration over the photon transverse momentum 
d^ qt using the well-known relation: c?^ Pi ^ — zfd^ qt However, the leptonic tensor depends not only 
on but also on cos 02- Therefore, factorization with the transformation to the integration over qt does 
not allow for calculating the angular dependence of the corresponding cross section. 



IV. FACTORIZATION PROCEDURE: LEADING ORDER as CORRECTIONS 

Now we proceed with the analysis of leading as corrections which are associated with the diagrams that involve 
gluon emission. It is known that large logarithms appearing in this case can be absorbed into the corresponding 
evolved fragmentation functions. We will pay a special attention to the terms with mass singularities which are 
extracted from the diagrams with the emission of real gluons. The corresponding Feynman diagrams with 0^5 
corrections are depicted in Fig. [21 The domain of integration over the loop momentum in each diagram is in 
the region where the considered parton is collinear to either Pi or P2 momentum directions. For the sake of 
definiteness, we assume that the hadron with momentum P2 belongs to the quark jet, while the hadron with 
the momentum Pi - to the anti-quark jet. The opposite situation is trivially obtained by the interchange of the 
labels. 



IVa. Evolution of the quark fragmentation function 

To begin with, let us consider evolution of the quark fragmentation function D{z). At the Born level, the 
quark which decays into the detected hadron with the momentum P2 and a bunch of undetected hadrons can 
emit a real gluon before and after interaction with the virtual photon. The diagram depicted in Fig. [5] (a) 
corresponds to the interaction of quark with the virtual photon before the emission of the real gluon. Choosing 
Feynman gauge for the gluon fields, we write the hadronic tensor corresponding to this diagram in the form 
resembling the Born diagram tensor: 

wW^''(FigIla)) = g'Cp J d% d^pS^^^ {ki + p - q) tr [7, e(fci) 7^ ^^(p)] • (39) 

However, instead of the quark correlator Q{k2), we have here the modified correlator fl{p) defined as 

^^P^ = ^^-l 0^([^^-^2ne(fc2)7a^. (40) 
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In principle, factorization of (|40p can be implemented in a similar way. However, in difference with the Born 
diagram, the modified tensor (|40)l in ((39)) is not completely soft because of the p-dependence. 

To factorize ([39]), let us first express the parton momenta in terms of the corresponding fractions of the hadron 
momenta by means of the integral representation of unity. The definitions of the fractions for parton momenta 
ki and k2 are the same as in ([T9| whereas for the loop momentum p we write 



d'^P-TT^Hp-n-l/y'). 



iv') 



(41) 



The pure soft part of (|1(I)) is associated with the quark correlator function Q{k2). At the same time, the quark 
propagators and the (5-function emerging from the imaginary part of the gluon propagator have to be associated 
with the hard part of (|40)) . Therefore, after expanding the (5-function in (|40|) about the momentum k2 around 
the direction defined by the hadron momentum P2, we get for the correlator: 



r!(p) 



■la 



S P 



2p-P2 



e(4)7a 



(42) 



From this tensor, which is the same as in the Born diagram, we single out the spin-independent quark fragmen- 
tation function (see eqs. (pS)) and (PS|) ). Inserting the tensor (|^^ into the hadronic tensor ([M)) . we get: 



>V(-)-(FiglIa)) = / 



dz[ f dy' 



dp^ 



i4r 



iy') 



p; 



y 



,4 ( vTTJxi [T"' (^1 ) 7m P 7a (^2 ) 7a P 




(43) 



-i/Q 



where we decomposed again the four-dimensional (5-function around the corresponding hadron directions, 



<5(4) (fci+p-g) 



f_2_ 

v' 



6^' 



Pi 



(44) 



and calculated the integral over dp with the (5-function coming from the imaginary part of the gluon propagator 
which fixes the minus component of the loop momentum: 



S{[p^k2]') 



1 



2[p+ P+lz',] 



Pa 



2[p+ - P+l-4] 



(45) 



We remind that p'^ and P2 are large vectors. Therefore, as one can see from (j45p . the minus component is thus 



suppressed as 1/Q and can be discarded in (|43)) . 

In (|^5)) . we still have integration over the plus and perpendicular components of the loop momentum. Just 
like in Fig. [5] (a) , the quark with momentum p emits a real gluon and transforms into the quark with momentum 
k2- This means that 



P+^^kt, dp+ = ^d^ 



(46) 



Using (|46|) and calculating the corresponding integral and the trace in ([39]) . we get the following contribution 
of the diagram [2] (a) to the hadronic tensor: 



wW^''(FigIla)) ^-2g^,^Ci.5(^) (Pi^) [ziD{z,)] 



'^"''■D{z'2) 



Pi 



(47) 



As in the case of the Born diagram, one can show that the 5-function 5(£^ — z'^j Z2) in (I47p should be associated 
with the imaginary part of the hard effective propagator, see [3] (a). Note that similar arguments are valid for 
other diagrams needed for the study of the quark fragmentation function evolution. 
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Integration over the two-component loop momentum should be implemented with the lower limits defined 
by the infrared cut-ofF and the upper limit of . Therefore, the factor ln((5^/A^) which appears after this 
integration reflects the coUinear singularity. 

Turning back to the averaged hadronic tensor, we thus derive 



W;,t^' '(FigUa)) = I d'Pii_ Wlt^ = -2 .9^ ^ Cf In 



as 



91 

A2 



Z2 



(48) 



The tensor (|48l) has a completely factorized form and corresponds to[3](a) which plays a role of a ladder diagram. 

The diagram [2] (b) does not contribute to the terms containing the mass singularity. Indeed, the hadronic 
tensor corresponding to such a diagram is given by 



>VW^'?(FigIlb)) = g'CF J d^hd''pS(''Hh+p-q)tr[Q{h)n^,{ki,p)] 



where 



..r _ , j3 + fci - fc2 



7m9(k2) -/u 7Z—r, TT^ 7a- 



(49) 



(50) 



In the same manner as before, one can see that, after integration over dp^ with the (5-function originating 
from the gluon propagator, the denominator of (|49p does not contain the term with that produces a mass 
singularity. As a result, this contribution can be discarded. It is necessary to note that the diagram [5] (b) 
effectively corresponds to the diagram Fig. [3] (b) with the self-energy insertion into the quark propagator. 
Two diagrams presented on [2] (c) and (d) contribute to the hadronic tensor in a following way: 



W^i^'^FigHc + d)) = 252Cj. J d^kid^pS^^^ (ki+p-q) tr[7,e(fci)r!^(fci,p)] 



where 



(27r 



^'^([^-^2] )7r-T— T-W 



[p + ki~k2Y p-^ 



(51) 



(52) 



In ([5T|). the denominator contains the necessary power of p^ and we have to keep only the zcroth order of p^ 
in the trace. Then, following the scheme outlined earlier for other diagrams, we derive the following expressions 
for the hadronic tensor 



^(^^■^(FiglStc) + (d)) = -25^, ^C^^ 5(2) (Pi^) [zi5(zi)] 
and for the averaged hadronic tensor: 



dz'2 



'-2 r,/ /N 24/^2 

(4)2 ^^^^^ 1 - Z'JZ2 



Pi 



(-L).9 



(FigIlc) + (d)) = -2.9^,^C7^.1n 



A2 



[ziD{zi) 



dz'2 

(4) 



2 4/z2 



4/^2 



(53) 



(54) 



Thus, we have all the ingredients for the derivation of the factorized hadronic tensor including the as corrections 
and mass singularities. Summing of eqs (I47p and (j53p together with the contributions related to the virtual 
gluon emission gives 



Wit^^' = -2g^^^CFln[:^][zMz,)] 



± as 



dz'o 



7T2^(4) 



(4? 



1 + (4/^2)- 

1 - Z^/Z2 



(55) 



The factorization scale fip can be introduced by the standard decomposition: ln(Q2/A2) — ln{Q'^/iij^) + 
ln(/x|,/A2) where the first term should be combined with the hard part of the corresponding hadronic ten- 
sor whereas the second one with the soft part. If we choose — Q^, the sum of ([55]) with the contribution of 
the Born diagram (I32|) leads to the following substitution for the quark fragmentation function: 



1 + 2/1 



y2 



(56) 
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As a result, the quark fragmentation function acquires Q^-dependence and satisfies the DGLAP evolution 
equation: 



^J^.I^Jlni..„,, v„,„,.||c,(i±^) ,57) 



IVb. Evolution of the anti-quark fragmentation function 



We will now consider the anti-quark fragmentation functions. The anti-quark sector in the fragmentation can 
be studied in similarly to the quark fragmentation function case. However, there are some minor differences in 
comparison to the quark fragmentation function evolution. First of all, the diagram [2] (b) will now play a role 
of the ladder diagram. Let us write down the hadronic tensor corresponding to this diagram: 



W^i'''(Figl2rb)) = g^Cp J d^k2d''m6^''Hk2+m-q)ti[^^e{k2)j,n{m)] 



where 



n{m) 



■la 



-^^([m-fci]2)e(fci)7„-2 



(58) 



(59) 



As in the case of the quark sector, we first introduce the definition of the corresponding parton fractions (see 
eqs. and (I4ip ). Then, the corresponding (5-functions have to be decomposed in the directions defined by 
the hadron momenta. We obtain 



J^'') {k2+m- q) 



y' 



V' 



(60) 



for the four-dimensional (5-function which is responsible for the momentum conservation at the local vertex: 



1 



2[m- - P^/z[] 



5 ni^ + 



2\m- 



(61) 



for the (5-function in the tensor ([59| . Integration over dm^ with J- function ((6T|) fixes the plus component of the 
loop momentum to be a small variable. As result, m"*" can be neglected in the corresponding expressions. 

Omitting the details, since all stages of calculations are exactly the same as for the quark sector. Below, we 
write the expression for the averaged hadronic tensor corresponding to Fig. [2] (b) : 



w!,t'''(FiglIb)) 



as 



^2 



dz'-^ 



5(4) 1-^ 



(62) 



The next non-zero contribution comes from the diagrams presented on Fig. [2](c) and (d). These diagrams give 
us the following expression: 



W(i)'«"(Figl2i;c-hd)) = 2g'CF J d^k2 d''mS^^Hk2 + m ~ q) tr [Oik^) -f,n^{k2,m)] 



where 



ri^(fc2, m) = — 7q 



d'^ki 



<5([m-A:i]2)e(fci)7^ 



rfi — ki + k2 



■7a- 



(63) 



(64) 



(27r)3 ' "''^'"^ (to- fci -hfc2)2 
Calculating the tensor (|63p in a way similar to the derivation of (|5ip . we obtain the averaged hadronic tensor: 



wit^'^"(Figlirc) + (d)) = -2g^, ^Cpln(^ 



Z2 



(65) 
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Diagram [2] (a) in this case does not contribute to the evolution of the anti-quark fragmentation function due to 
the same arguments as we presented for the quark sector discussing the diagram [2] (b) . 

Thus, combining aU the diagrams and adding the contributions from the virtual giuon emissions, we get: 

Again, as for the quark case, the summation of the hadronic tensor (|66p with the Born hadronic tensor modifies 
the anti-quark fragmentation function: 

zi 

As a result of it, the anti-quark fragmentation function becomes Q^-dependent and obeys the DGLAP evolution 
equation (see, ([57])). Besides, the introduction of the hard effective propagators as in the quark case leads to 
the factorized Feynman diagrams which are completely analogous to the diagrams on Fig. [3l 



1 + 

1 - z'Jzi 



I 




[^(^2)1 






^2 



(66) 



V. /CT-DEPENDENT FUNCTIONS: BORN APPROXIMATION 

In the preceding sections, we focused on the case of the spin-independent integrated fragmentation functions 
which in a way was a test of our approach. Since we plan to extend our approach to fcy-unintegrated functions (in 
particular, Collins function) we discuss the Born approximation of e+e^-annihilation, involving both the Collins 
fragmentation function and the spin-dependent fragmentation function. Namely, we consider the production of 
transverse polarized quark-antiquark pair associated with (T^^-structures in the corresponding matrix elements, 
see below. 

More exactly, the nonperturbative blob (see, for example. Fig. [1]), related to the detected baryon, whose 
transverse polarization is correlated to the transverse polarization of the quark can be described by the spin- 
dependent fragmentation function. At the same time, in the other blob, the transition of the transverse polarized 
antiquark with the intrinsic transverse momentum into the unpolarized hadron is described by the Collins frag- 
mentation function. In latter case, the mentioned transition is related to the azimuthal asymmetry distribution 
of hadrons. Note that the symmetric case, with two similar fragmentation functions, has been well-studied from 
both the theoretical and experimental points of view Bin, is |21| . The asymmetric situation requires the 



simultaneous measurement of the meson azimuthal asymmetry in one jet and baryon (say. A) polarization in 
another one, Such measurements may be performed at BELLE [2^ . 

We use the co-ordinate (or the impact parameter) representation where the explicit definition of the transverse 
momentum is not necessary. In this case, the corresponding hadronic tensor takes the form: 



7Mei(07.©2(0 ■ (68) 



Here, the spin-dependent fragmentation function has been defined in the co-ordinate space as 

1 

QiiO j ^e-(^-«)/^^i/iT(4), (69) 



where S denotes the hadron transverse polarion, and the Collins function's analog in the coordinate space is 



given [15| by 

1 

0^(0=^ ^Ma^0P,-eJ^e-^^^^<y<Ht{z[), Ht{z[) ^ j dk?, H^{z[,kll (70) 


where ^ is the (transverse) position in co-ordinate space; while M, introduced due to the dimensional analysis, 
implies a parameter of the order of jet mass being the only dimensionful parameter in the soft part. 
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The position in the coordinate space ^ (see, (j70p ) yields the derivative over Pi in momentum space. Therefore, 
inserting §9^ and ^ in we get: 



dz'i f ( d 



Pi 



where 



Siq- - ^) }S{q+ - ^)S^^\-L)H^iz[)HiTiz'2)Tf..{Pi,P2,ST), (71) 



(72) 



T^.{Pi,P2,St) = tr [7^(iM(7+-Pf )7.(a-T75P2+^T)] • 
Using the kinematics defined above, we derive 



dz[ {z[Ht{z[))5^^\z[-zi) 



^i/lT(4)'5(4-^2) 



Then, calculating the averaged hadronic tensor, one has the following: 



Hit{z2) 



Z2 



(73) 



(74) 



Note that the analogous asymmetric combination corresponding to SIDIS process related to ours by crossing 
was also considered in [15'| . In that case, the incoming quark in SIDIS is described by the transversity distribution 
function: 



(75) 



where S'^ is the target polarization. Using ([70]) and ([75)) . the corresponding hadronic tensor takes the form: 



(27r)' 



(76) 



To study the kx (or P2±) distributions in SIDIS and the related asymmetries, we should consider the weighted 
hadronic tensor which projects out the corresponding moment of the Collins function: 



A„W 



d^P25{Pi){P2-ni_)5 



Pi -Pi 
Pl-q 



z Aw'^^ 



(77) 



where n± is the unit transverse 4— vector {n± ■ Pi = n± - q = Q, n\ = — 1) which defines the transverse direction. 
Using (|70p. ([75|) and ([77]) . one can see that the derivative 9" = djdP^ in 



A„>V" = iM / d*F2 {P2 ■ nx_) (5(P2 ) ^ 



Pi ■P2 
Pi-g 



dxdz'd°'5{xPi +q- P2/ z')h{x){z' I{z'))iiYi5PiS-f''[P2-ic.]Y\ 



(78) 



should act only on (P2 • rij^), so that the result [iBj is equal to the standard expression for the contribution 
of Collins function, except that the role of the direction of intrinsic transverse momentum is played by the 
auxiliary transverse vector n^: 



(79) 



This substitution does not change the azimuthal dependence, as the weighted integration corresponds to az- 
imuthal average: 



{d(j{(f)h)cos{(j)h ~ (j)n)) = cos(/)„(dcr((^ft)cos(0,j)) +sin0„(dcr(0,,)sin(0,,)). 



(80) 
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As a result the azimuthal dependence of the cross-section is transferred to the dependence on the angle and 
I{z) corresponds to the moment of the Colhns function: 

I{z)^ j dk?r^H^{z,kl). (81) 

Thus, to describe the Colhns effect we suggest the /c^-dependent fragmentation function to be written in 
the co-ordinate space, where no specification of intrinsic Ut is required. Therefore, the calculation of radiation 
corrections and evolution may be performed in the same way as in this paper. 

Moreover, this approach may be applied [l5| for the description of higher weighted moments in fcy and pr- 
Say, choosing the term of order in the expansion of unpolarized fragmentation function in the coordinate 
space and taking the pT moment of the hadronic tensor weighted with p^ corresponds to the account of the 
width of the fcr-dependent unpolarized fragmentation function. The exponential shape of the latter, in turn. 



corresponds to partial resummation 15| of the infinite series of higher twists, analogous to that leading to 
the appearance of non-local vacuum condensates, when vacuum rather than hadronic averages are considered. 
The expansion in coordinate space accompanied by taking the respective weighted moments of cross-sections 
provides a complementary definition of observables, accounting for the shape of /cr-dependent fragmentation 
(and distribution) functions. 



VI. CONCLUSIONS 



We described a method which allows us to prove a factorization of the process with two fragmentation 
functions. We would like to point out that the presented method can be applied for any two-current process. The 
difficulties of factorization for such kind of processes emerge for the case when the kinematical transversalities 
inside the hadron are rather small. It leads to the problem in the definition of what is the hard subprocess for 
the process. 

Following the idea of the paper IJ] , we showed that the corresponding (5-functions in the hadronic tensors 
should be treated as the hard parts. It is based on the observation that these (5-functions can be associated 
with the imaginary parts of the effective propagators related to the well-defined hard subprocess. As a result, 
we finally have the completely factorized expression for the hadronic tensor with the evolved fragmentation 
functions. 

In this paper, the proposed method has been tested in the simplest case when the differential cross section 
of e+e^ annihilation is related to the spin-independent integrated fragmentation functions. We also extended 
our approach to the study of the spin-dependent structures and fcT-dependent fragmentation functions (Collins 
function and transversity fragmentation function). This will allow us to perform the leading order QCD fits of 
relevant experimental data [19| . 
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Figure 1; Non-factorized (a) and factorized (6) Born diagrams 




Figure 3: Factorized leading order as diagrams. 



